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SECOND-ORDER SMALL-DISTURBANCE SOLUTIONS FOR 
HYPERSONIC FLOW OVER POWER-LAW BODIES 


James C. Townsend 
Langley Research Center 

SUMMARY 

Similarity solutions have been fovmd which give the adiabatic flow of an ideal gas 
about two-dimensional and axisymmetric power-law bodies at infinite Mach number to 
order-6 , where 6 is a body slenderness parameter. These solutions were obtained by 
a perturbation method wherein the flow variables are expressed as series of similarity 
functions and the axial variations in the flow equations separate out. The resulting simi- 
larity equations were numerically integrated from the shock wave, where the boundary 
conditions are known from the Rankine -Hugoniot relations, toward the body. The 
order-6 solutions are limited to body power-law exponents greater than 0.85 for two- 
dimensional flow and 0.75 for axisymmetric flow to avoid a singularity in the shock- 
displacement calculation. The solutions, which are independent of the slenderness 
parameter 6 and, thus, are universal functions, are presented in tabular and graphic 
form. Because of the entropy layer induced by the nose bluntness (for power-law bodies 
other than cones or wedges) the velocity and density functions are not valid at the body 
surface, but the pressure function is unaffected. 

In comparisons with the exact solutions for inviscid flow over wedges and cones, 

O 

the order-6 similarity results give excellent agreement for 6 values up to 0.4, corre- 
sponding to wedge or cone angles up to about 20°. In comparisons with experimental 
shock-wave shapes and surf ace -pres sure distributions for 3/4-power axisymmetric 
bodies, the order-6^ similarity solutions give good results, considering that Mach number 
and boundary -layer displacement effects are not included in the theory. For body fine- 
ness ratios near 2, the effects of the order -6 terms are significant only very near the 
body nose, whereas for a fineness ratio near 1 the order-6^ terms have a large effect 
over almost the entire body. 
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Although the order-6 similarity solutions are developed for infinite Mach number, 
the derivation shows that they are compatible with shock-strength perturbation solutions, 
which introduce Mach number effects. 



INTRODUCTION 


A great deal of research has gone into finding solutions to the hypersonic small- 
dist;irbance form of the inviscid adiabatic -flow equations. One area that has received 
particular attention is that of similarity solutions for power-law profile bodies. (See 
refs. 1 to 6.) Since these studies have all begun with the small -disturbance equations, 
obtained by neglecting terms containing the body thickness ratio, their solutions also 
neglect the effects of body thickness. The purpose of the present study was to obtain 
second-order similarity solutions for power-law bodies by retaining all second-order 
terms in the body -thickness ratio. (As will be seen, there are no first-order terms in 
the body thickness and, thus, no first-order solutions in the thickness parameter.) 

The similarity solutions referred to herein are solutions for self -similar flows, 
that is, flows for which the flow field between the shock wave and the body can be 
expressed in terms of functions which (in suitable coordinates) are independent of one of 
the coordinate directions. Inviscid axisymmetric supersonic flow over a cone with an 
attached shock wave is a classical example of a self-similar flow and represents a par- 
ticular case of the similar solutions discussed herein. For the cone, the flow -field prop- 
erties (e.g., the pressure, the density, and the velocity components) are themselves con- 
stant along rays from the cone vertex. For the other power-law bodies, the flow -field 
properties are not constant themselves, but similarity functions describing these prop- 
erties are constant (to the order of the solution) along curved power -law paths from the 
nose of the body. The similar-solution approach to solving the flow equations is valuable 
because it allows a reduction in the number of independent variables in the problem. In 
particular, for hypersonic flow about power-law bodies, the similarity approach reduces 
a system of partial differential equations to a system of ordinary differential equations. 

Kubota (ref. 3) obtained numerical solutions for the self-similar flow fields about 
bodies having power-law profiles in the limits of infinite Mach number and very slender 
bodies (the "zero-order" case). He also applied a perturbation in a shock strength 
parameter and numerically obtained first-order similar solutions for moderately strong 
shock waves. Mirels (ref. 4) computed additional and more accurate numerical results 
for the zero-order and moderately strong shock-wave cases. He also derived approxi- 
mate analytical solutions for these cases. Hayes and Probstein (ref. 5) described the 
general development of similarity solutions and emphasized the contributions of foreign 
authors. Mirels (ref. 6) treated all of the important developments in the use of h 5 rper- 
sonic small -disturbance theory to obtain solutions for power -law bodies in a unified way, 
and added an analysis of perturbed power-law body shapes. 

Only a few experimental investigations of the flow field over power-law bodies have 
been made. Kubota (ref. 3) compared his theoretical results to measurements of 
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siirf ace-pressure distributions and shock -wave shapes for 2/3- and 3/4-power bodies 
and obtained good agreement for the more slender bodies. Peckham (ref. 7) measured 
pressure distributions and shock-wave shapes for a series of power-law bodies, some of 
which fall in the similar -solution rai^e. Freeman, Cash, and Bedder (ref. 8) and 
Beavers (ref. 9) also presented detailed shock-shape data for a series of power-law 
bodies and registered some disagreement with Kubota’s results. Spencer and Fox 
(ref. 10) presented aerodynamic drag and other data for several power-law bodies over 
a wide Mach number range. Ashby (ref. 11) presented aerodynamic data for a similar 
series of bodies over a range of Reynolds numbers at a Mach number of 6, and Ashby 
and Harris (ref. 12) used a method of characteristics and boundary-layer computer pro- 
grams to show the important effect of boundary-layer transition on the total drag of those 
bodies. 

Townsend (ref. 13) applied the zero-order solution of Kubota and Mirels, with their 
shock-strength parameter perturbation and with a boundary-layer displacement correc- 
tion, to the problem of estimating the forces and moments on a half-axisymmetric body 
under a thin, flat wing. In order to study a range of configurations at a moderately hyper- 
sonic Mach number, Townsend applied his method to configurations which are marginally 
slender (i.e., to configurations for which the errors arising from body thickness are 
small but not negligible). This type of application points up two reasons for seeking 
solutions which include the effect's of the second-order terms for body slenderness in the 
flow equations; (1) to assess the error caused by making the small -disturbance assump- 
tion, and (2) to improve the accuracy of calculations for marginally slender bodies. 

In the present study the second-order similarity solutions were obtained by a per- 
turbation method. This method used expansions of the variables in terms of the body 
slenderness parameter to obtain higher order solutions as perturbations from the known 
zero-order solution. This approach is very similar to that of Kubota and Mirels in 
obtaining their first-order solutions in terms of a small strong-shock parameter. A 
more detailed presentation of the theoretical development is given in reference 14. 

SYMBOLS 

^l’^2 shock-wave displacement constants (eq. (6)) 

p - p 

Cn pressure coefficient, ^ 

^ q 

^OO 

D length constant in shock-shape correlation (fig. 11) 
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body fineness ratio, 2 

Maximum thickness 


similarity pressure function 

body length 

power -law exponent 

free -stream Mach number 

pressure 

free-stream dynamic pressure, - p 

2 

lateral coordinate 

lateral coordinate of shock wave 

temperature 

longitudinal -velocity perturbation 
free-stream velocity 
lateral velocity 
velocity 

longitudinal coordinate 



00 CO 


alternate shock-shape parameter, 
ratio of specific heats 


2(1 - m) 
(1 + CT)m 


Rn(r) 

body slenderness parameter, — ^ — 

I 

2 

shock-wave strength parameter, (SM^) " 



v 


similarity lateral coordinate 


0^ body surface angle 

0g shock -wave angle 

0Q similarity stream function 

V similarity longitudinal -velocity function 

^ similarity longitudinal coordinate 

p density 

cr constant with value 0 for planar flow and 1 for axisymmetric flow 

(p similarity lateral -velocity function 

\jy similarity density function 

Subscripts: 

a part of decomposed function multiplied by a .2 (eqs. (22)) 

b body surface 

c part of decomposed function independent of R 2 (eqs. (22)) 

s shock wave 

t total 

0 zero order 

1 order -€ 

2 order-6^ 


OO 


free stream 


A prime denotes the derivative of a function of one variable. A bar over a variable 
denotes that it has physical dimensions. Normalized variables (without a bar) are given 
in equations (1). 


THEORY 


Transformation of Basic Flow Equations 


Normalization. - The initial treatment of the basic flow equations follows that of 
Kubota (ref. 3) (also covered by Mirels in ref. 6), except that no terms are dropped. 
They showed that for slender bodies in hypersonic flow the variables can be normalized 
by using the expressions 




x = S 


r =4 

51 


P = 


00 ^ 00 ^ 




U - u. 


u = 


s2- 

6 u 


V = 


6u 


OO J 


( 1 ) 


The symbol 6 is a body slenderness parameter (to be discussed later) introduced so as 
to make the dimensionless variables of order unity. These variables are substituted 
into the complete inviscid adiabatic -flow equations to obtain the following normalized 
flow equations: 


Continuity: 


dx dx dr r 


Longitudinal momentum: 


62uM*|i + v|S + i|E=0 

8X dx dr p ax 


Lateral momentum: 


Energy: 


dx dx dr p dr 


^ dx dx dr)[py) 


J 


( 2 ) 
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The constant a in the continuity equation has the value 0 for planar flow (Cartesian 
coordinates) or the value 1 for axisymmetric flow (cylindrical coordinates). Note that 
each of these equations contains a leading term in 6^. If the body were sufficiently 
slender, the order-6^ terms could be dropped, thus leaving the hypersonic small- 
disturbance equations used by other workers. For the present study, however, the equa- 
tions are retained in the complete form. 

Similarity variables .- The next step is to put the flow variables into similarity 
forms. Following Kubota's method, these forms will be selected so as to agree with the 
flow through the oblique bow shock wave. The normalized flow variables just behind an 
oblique shock are (ref. 15) 


» on 


I 2yMj sin^ds - (y - 1) 


y + 1 




(y + 1)M^^ sin^e^ 

(y - 1)M 2 sin20„ + 2 




■-2( 


2(M„^ sin^Sg - l) 


(y+.i)Mj J 


Vc = - 
® 6 


2{mJ sinh^ - l) 


L (y + 


cot 0, 


(3) 


If the shock -wave shape is given by R(x), its slope is R' (x) s ^ = 6^ = 6R*. The 

dx dx 

shock -wave angle 0g is related to the slope by tan 0g = 6R', from which 
. 2 6^R'2 

sin 0g = g — JJ-. Putting these results into the oblique -shock relations gives the 

1 + 

pressure (for R' of order unity); 


2 

r r'2 

- y - 1 g 

y + 1 

1 + 62r'2 

2y 


2 

y + 1 


R’ 


.2 


— ^ 62r'‘* 
y + 1 


2^i + o<64) 

y + 1 y 


(4a) 
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Similarly, the other flow variables at the shock are 


Ps 


y + 1 _ 2(y + 1) e_ 

y - 1 (y - 1)2 R»2 


+ 0 ( 64 ) 
» 


(4b) 


u 


s 


V 


s 


R’2 + _A_ e + 0(64) 

y + 1 y + 1 y + 1 

_1_. R* _ _2_ 62r, 3 _ _2_ _e_ ^ 
y + 1 y+1 y+lR' 


(4c) 

(4d) 


where e s is a shock-strength parameter; and as e 

62m 2 

oo 

limiting value for shock -wave strength. 



y + 1 
y - 1’ 


the 


At this point, in order to get an expression for the shock -wave shape, consideration 
is narrowed to flows about power -law bodies. Under the hypersonic small -disturbance 
assumptions, a power-law body ^r^ ~ produces a power-law shock wave (r ~ x”^) 

for — - — < m S 1. (See ref. 2.) Specifically, for 6^—0 and e — 0, a body 
3 + or 

rjj = produces the zero-order shock shape Rq = x”^. (See figs. 1(a) and 1(b).) 

Note that for m = 1 the body is a wedge (for a = 0) or a cone (for ct = 1), both of which 
are known to have straight shock waves and therefore satisfy the aforementioned relations. 


The expression for the zero-order shock-wave shape serves to define the slender- 
ness parameter 6. That is, Rq(x) = 6f 

6 = — = — . Thus, 6 is the tangent of an angle defined by the shock -wave position 

(fig. 1(a)) and is, in fact, a "mean shock-wave angle" parameter. However, the shock 
lies near the body so that the shock-wave angle and body slenderness are closely related. 

To aid in the separation of variables, a shock -oriented coordinate system is now 
introduced (fig. 1(c)). This system has 



evaluated at x = f gives the relation 


^ = x 



( 5 ) 


so that r = tjRq = r]x^ = 


The body surface is then 


^b = where 
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Lateral coordinate, -n Lateral coordinate, r Lateral distance. 


S/h 



Longitudinal distance, x 

(a) Physical coordinate system (shown with 6 = 0.4). 



Shock, 17 = 1 


Body. 97 =77^ 





Longitudinal coordinate, ^ 

(c) Similarity coordinate system. 


Figure 1.- Power -law body and zero-order shock in physical 
and transformed coordinate systems (shown with m = 0.85 
and a = 0). 
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The shock -wave shape to be used in equations (4) is the zero -order shock with 
Kubota's shock-strength perturbation (ref. 3) and a separate perturbation for the body 
slenderness. It is taken to be 


R(^) = + 6^a2m2^"2(l-“) 


( 6 ) 


where the exponents are of the form required for similarity and the constants a^^ 
and a2 are to be determined as part of the solution. 

Substituting the shock -wave slope (derivative of eq. (6)) into the expanded oblique - 
shock relations (eqs. (4)) and ordering the terms by powers of 6 and e suggest 
expressions for the flow variables in the following form (neglecting terms of order 6“^ 
and of order e2); 


= FQ(77)m24-2(l-m) + eF^{n)m^ + 

pU,V) = ^{^oiv) + 






(7) 


J 

Here, functions of 77 = replace the constant coefficients of the shock relations. 

Transformed flow equations .- Equations (7) are substituted into the normalized 
flow equations (2), noting that 


JL=1L~ ^ — m77 9 

dX dX dx dr] ^ dr) 

} 

d_d^d^dr]d_l d 
dr dr 94 dr dr] dfj 


( 8 ) 
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Then, the flow equations (2) become (away from the nose ^ = 0): 
Continuity: 

[*^0 + f ‘^ o ) - - ^ o ) ^ {^ 0^1 ^ [^0 + f >^ 0 ) 

+ [^ 0 ' + -{v - 0o)’^l} 

+ + (’^o' + f "^o] ^2 ^ ^ *^0 ~ ^(Sn^’^ 2 " “ *^0) ’^2' 

- ^ ^0*^0' - [^ 1^0' + ^ = ° 


Longitudinal momentum: 


(i-i'oH' 


—){''o*° 

m y \ u 




0/J 


„3j.-3+2m 
m ^ 


+ e< 




+ 6^^ (77 - 0o) +Jl^+ - '^0'^2 - ;^[^^0’ + 


^o (5 V - 2(4^2) 0(64) . 0 


(9b) 
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Lateral momentum: 



These are ordinary first-order differential equations that are linear in the derivatives of 

the functions defining the pressure, density, and velocity fields. It is noteworthy that 

although 6 appears to the first power in the normalization of variables (eqs. (1)), only 

the even powers of 6 appear in the final form of the flow equations. Thus, although the 

solutions to be found are of second order in the body slenderness parameter 6, they 

2 

could be considered of first order in 6 . To avoid any ambiguity they will generally be 

O 

referred to as order-6 solutions. 
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Similarity equations.- Since the general similarity solutions being sought do not 
depend on the particular values of 6 or e, each of the three major terms in each of 
the four conservation equations (9) must be separately equal to zero. Therefore, the 
terms can be separated into 12 equations in the 12 imknown fimctions Fq, x^q, vq, 

(Pq, Fj^, (py^, Fg, xf/2> ^2‘ these equations form the zero- 

order and order-e systems of equations foimd by Kubota (ref. 3). 

The zero-order equations are as follows; 


Continuity: ^ 

f *^o) ■ " H) ° 

Longitudinal momentum: 

Lateral momentum: ) 

(” - ■'■o) V - ^ *0 " “ 

Energy: 



( 10 ) 


These equations are the same as the case first studied by Kubota, except that he omitted 
the longitudinal -momentum equation, which is uncoupled from the others. These equa- 
tions contain only one parameter, the power-law exponent m. Thus, for two-dimensional 
flow (ct = 0) or axisymmetric flow (cr = 1) of a given gas, the similar solutions Fq(t 7 ), 
xpQ{rj), Vq{t]), and <^q(tj) each form two families of "universal functions" depending only 
on the power law of the body. 
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The order -e system of equations is given as follows: 


Continuity: 




Longitudinal momentum: 




Lateral momentum: 


(•) - ^ ^ '^1 - [^ 0 ' + (Hr)]*i = “ 


(lib) 


(11c) 


Energy: 


{1 - 0 o ) 




*0 "0 


^0 


2 n -:^"i, 

0 y 




I *'0 ^oP 


(lid) 


These equations are the same as Kubota’s (ref. 3) first-order perturbation for shock- 
wave strength, except that (again) he omitted the longitudinal -momentum equation since 
it is uncoupled from the rest. They can be solved numerically by using the results 
of the zero-order solutions. The resulting similarity functions Fj^(t;), v-^{r]), 

and 0j(r;) are also universal, like the zero-order fimctions, in that they depend only 
on m as a parameter. References 3,4, and 13 contain the results of numerical solutions 
of the zero-order and the order-e equations. 

The remaining four equations form the order-6^ system, which has not been deter- 
mined previously, and are given as follows: 
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Contimiity; 


^ 0^2 + {^0 + ? ^o) <^2 - - ^ 0)^2 ^ [^ O ’ + ? - KHr )]’^2 

- - [^ V + ° 


(12a) 


Longitudinal momentum: 


I” - *o) •' 2 *^^* HHr) (■'2 * ^) ■ 

- %' 4>2 + *'o[^V + ° 


(12b) 


Lateral momentum: 


{■J - t'o)02 - IT * A -^2 - k' - 3{kr)l'^2 *ko' (kr)^0 


1^0 = 0 


(12c) 


Energy: 

{V - cPo) 


^ 2 ' ^2*\ ( ^0^2 ^o’^2' 


*0 I'd 




0 


'' \*o *'oy 


■('¥-^>.- N - K-^l 


Vq = 0 


(12d) 


Except for additional terms corresponding to the order-6^ terms retained in the normal 
ized flow equations (2), these eqviations are very similar to the order -e equations; many 
of the coefficients are the same, and the only body shape parameter that appears is m. 
The similarity functions F 2 (t 7 ), ^2^^’ which form the solutions to 

these equations, will, therefore, be families of universal functions in the same sense as 



the other solutions are. Furthermore, just as the order -e equations are independent of 
the body slenderness perturbation (in these equations are independent of the shock- 
wave strei^th perturbation (in e). Thus, application of these equations to determine the 
body slenderness perturbation of the ’zero -order small -disturbance similar solutions 
neither requires nor excludes application of the equations for the shock-strength pertur- 
bation at the same time. 

Since the order-e solutions have been fovmd previously and are not needed to get 
the order-6^ solutions, they will not be considered further. (Ref. 14 shows that the 
order-6^ and order-e boimdary conditions are also mutually independent.) All subse- 
quent development will assume e « 6^ so that e = 0(6^) « 1; all terms of order-6^ 
or smaller will be neglected. 


Boundary Conditions 

Shock wave .- The boimdary conditions at the shock wave are determined by the 
oblique-shock relations (eqs. (3)). By usii^ the expression for R(|) (eq. (6)), these 
equations in their expanded form (eqs. (4)) become (omitting the order-e terms): 


® y + 1 y + 1 


^ C”ni +0(64) 


S y - 1 


Ug m^^ 

y + 1 


2 ,^2t-2(l-m) _ .2 




) (13) 


4,-4(l-m) ^ 


2 ^^-(l-m)^,2^ 

“ y + 1 y + 1 


(2S^)a2 - ^ o(64) 


Comparing these equations term-by-term with equations (7) determines the boundary 
conditions for the similarity functions at the shock wave (ry = Tg): 



'^O^^S^='y+l 

■''oOa) = ^ 

^ 1 


(14a) 
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^2^Vs) = 0 

- l] 

^ 2 ‘”=> = 7^[(^)“2 ■ 3 

Note that the shock -wave displacement constant a .2 initially is unknown. It depends on 
the parameter m and is to be found in satisfying the boundary condition at the body sur- 
face as part of the solution of the flow equations. 

Body surface .- The boundary conditions at the body surface are determined by 
requiring zero mass flow through the surface. Thus, the mass flow normal to the sur- 
face (fig. 2) is given by 

PbVb cos sin = 0 (15) 




Figure 2.- Vector diagram of flow at the body surface. 
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Or, in the normalized variables, 


Vb = |(l + 6\)tan 0^, 

Now, 


( 16 ) 




Putting this into equation (16) along with the expressions for u and v in terms of the 
similarity functions from equations (7) gives the relation 

The resulting boundary conditions at the body surface are 


^oK) -% = ^ \ 

02(^b) - 0 ^ 


(17) 


SOLUTION OF EQUATIONS 

The similarity forms of the flow equations (10) and (12) with the boundary conditions 
at the shock wave (eqs. (14)) and at the body surface (eqs. (17)) are sufficient to determine 
completely an order-6 solution for the flow field. Although there is no general analytic 
form, solutions can be foimd by integrating numerically for each value of the power -law 
exponent m (with cr set equal to either 0 or 1). This section first puts the equations 
into the form required for the numerical integration; then, it develops methods for over- 
coming certain difficulties in applying the boimdary conditions; and, finally, it describes 
the numerical -integration method. 

General Scheme of Solution 

The general scheme for obtaining the solution is to begin at the shock tj = 77 where 
the boundary conditions are known, and to integrate the similarity functions numerically 
toward the body, which is known to be reached when the zero-order boundary condition is 
satisfied; that is, V =V]j when 0q(t7) = 77 . The derivatives of the similarity functions, 
used for the integration, are found by solving the flow equations for them algebraically. 
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Thus, from the zero-order equations (10), the derivatives are 



The derivatives of the order-5 
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functions, from equations (12), are 
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(19b) 



(^7 - 





+ 


(^Q + ^ >/^o) ^ (^ir) ^ 0*^0 " ^ {%^0 



(19c) 



Two major difficulties must be overcome in order to apply the scheme of inte- 
grating these equations from the shock to the body. One difficulty is the singularity at the 
body surface apparent from the fact that the denominators of some of the terms of these 
equations approach zero as the independent variable 77 approaches the surface value 7]^. 
(From eqs. (17), <Pq - V ^ as rj -* It is overcome for the zero-order functions 
by using an approximate analytic solution developed by Mirels (ref. 4). Mirels' results 
at the body surface can be expressed as 



j /(t? - <^o) 

\v J y - ^\ 'n ) 

^ 2(2y - 0) Fq\71 I 


*'o(’'b) ” - I’ibC” - I 

M %) ' “ 

2 O’ 

These relations are correct to the order of 6q , where 6q = ~ <Pq) is a 

stream -function similarity variable. Since 6 q 0 as 17 — 77 ^, they apply near the body 
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surface. The order-6 similarity functions F2, ^ 2 > ^2 calculated at 

the body by extrapolation. The extrapolation curve (chosen for convenience) is the cubic 
in 77 passing through three of the computed points of the function and having its second 
derivative equal to zero at the body surface. The points used to define the cubic were the 
computed point nearest the body surface and two previous points chosen so as to make 
the distance between points larger than the distance between the last point and the body. 

The second of the two major difficulties is associated with the fact that the problem 
is a two-point boimdary -value problem. This difficulty is manifest in the need to choose 
the correct value of the shock -wave displacement parameter a2 (eq. (6)) at the beginning 
of the integration in order to satisfy the order-6^ boundary condition at the body surface 
at the end of the integration. The measures taken to deal with this difficulty are described 
in the following section. 

Methods for Determining the Constant a£ 

Since the constant ag is initially vmdeter mined, the value of the order -6^ shock- 

wave similarity ordinate ?7g = 1 + is unknown and cannot be used to 

begin the integration toward the body surface. Also, using the shock ordinate would 

reintroduce the longitudinal distance which is undesirable. The use of the zero- 

order shock ordinate 77 = 1 as the starting point for the integration avoids these two 

problems but requires that the boundary conditions be transferred from 77 = 77 g, where 

they are known, to 77 = 1 . This transfer is made by using the Taylor series expansions 

of the similarity flow variables about the point 77 = 77 . Applying the Taylor series 

s 

expansion to equations ( 7 ) (with terms of order e neglected) and evaluating at 77 = 1 
2 

gives to order 6 : 

P(«,l) = iPo(’’s) + ®^[^ 2 {’>s) - 
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Comparing the results of evaluating these equations (by using eqs. (18) for the deriva- 
tives and eqs. (14) for the boundary conditions at 7]g) with the results of evaluating 
equations (7) at 77 = 1 directly gives the transferred boundary conditions at 77 * 1: 



Iteration method.- These transferred boundary conditions provide a definite 
starting position for the integration toward the body, but the constant ag must still be 
determined by one of two methods. The more obvious one is to guess the value of a. 2 , 
integrate toward the body (using the method given in the previous section to reach the 

p 

surface), test the order-6 boundary condition at the surface, and repeat using improved 
guesses until the surface boimdary condition is satisfied closely enough. Improved 
guesses for this iteration method were made by using the method of chords, a finite dif- 
ference approximation to the well-known Newton-Raphson method. 

Decomposition me thod.- The other method for determining a .2 takes advantage of 
the linearity of the equations in the order-6^ functions, which allows superposition of 
solutions. It was used by Kubota (ref. 3) and Mirels (ref. 4) in obtaining their results and 
is applied in a similar mminer here. Each of the order-6^ similarity functions is 
decomposed into a linear combination in the parameter a£; for example, 
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= 4'2^J~'n)a.2 + ^2 

J 

9 

These decomposed variables are then substituted into the order-6 flow equations so that 
the continuity equation (12a), for example, becomes 



Splitting each of the equations obtained in this way into two separate equations (by setting 
the term containing the othei* term each equal to zero) produces a system of 

equations in the subscript-c functions and a system in the subscript-a functions. The 
system in the subscript-c fvmctions is identical to the original system of equations (12). 
The system in the subscript-a functions is the same except that the inhomogeneous 
terms (i.e., the terms that do not contain an order-6^ function or its derivative) do not 
appear. These two systems of equations have different sets of bovmdary conditions. In 
order to obtain them, the bovmdary conditions at 77 = 1 (eqs. (21)) are decomposed by 
comparisons with equations (22), thus giving 


F2,a(l) = 

" (i 

y + 1 [V 

3 m - 2\ 
m / 

(^)l 

f2y - 

( y - 1 


^2 c^^) ^ 

y + 1 
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y + 1 
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02,c^l) = 


2 

y + 1 


Beginning at 77 = 1 with these boundary conditions, the decomposed system of equations 
(which is now 12 linear differential equations in 12 unknowns) is integrated toward the 
body. The boundary condition at the surface, expressed in terms of the surface values 
of the decomposed functions, is (from eqs. (17)) 

Solving for ag gives 

- (23) 

■f2.a(->b) 

This value of a -2 is now used to recombine the decomposed similarity functions by using 
relations such as equations (22). 

Once these functions have been computed for any value of the body power-law 
exponent m (with ct = 0 or 1), they can be used to calculate the complete flow field 
about any such body as long as it is slender enotigh that 6^ « 1 and the Mach number 

is large enough that 6 = — L_ « 1. The appendix shows how to apply the solutions to 
a given power-law body. 


Description of the Numerical Method 

Two separate computer programs were written to integrate equations (18) and (19); 
one uses the iterative method and the other uses the decomposition method for obtaining 
the value of a 2 ~ These computer programs use a standard integration subroutine 
employing the fourth -order Runge-Kutta formula supplemented by a Richardson's 
extrapolation. This subroutine halves or doubles the integration step size automatically 
in order to meet specified local truncation-error criteria. 

For the present computations the initial step size (in 77), which was also the chosen 
maximum step size, was 2~'^ (0.0078125). Generally, the step size decreased to less 
than 2"15 near the body. At each step, estimates of 77^ and Fq^? 7^^ were computed 
by use of equations (20). When both estimates agreed to within 1.0 x lO"^ on successive 
steps, the estimates were accepted as the actual values of 77^ and Fo(^b) 
values of the other functions at the body were computed by the cubic extrapolation. The 
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maximum extrapolation distance for m ^ 1 was 0.000174 in rj. The iteration for ag 
was considered to have converged when the order-6^ boundary condition at the surface 
(eqs. (17)) was satisfied to within 0.5 x 10“10. 

DISCUSSION OF RESULTS 
Zero -Order Functions 

The methods given in the previous section have been used to compute the zero- 

p 

order and order-6 similarity functions for a number of cases. (All are for y = 1.4, 
representing an ideal diatomic gas.) The figures presenting the functions were generated 
by plotting machines directly from the computed results. The slight waviness which may 
be noticeable in some of the figures is a result of this computer-aided plotting process; 
however, the curves at all points on the plots are accurate to within ±0.1 percent of the 
fvill -scale values. 

The zero-order similarity functions Fq, i//q, 0q, and Vq are shown for several 

values of the power-law exponent m in figure 3 for two-dimensional flow (a = 0) and in 
figure 4 for axially symmetric flow (ct = 1). Numerical values of the ftmctions are given 
in tables I and n. These functions agree with the same functions calculated by Kubota 
(ref. 3), Mirels (refs. 4 and 6), and Townsend (ref. 13). 

The pressure function Fq and the lateral-velocity function <pQ are seen to be 
smooth and well behaved from the zero-order shock location {ri = 1) to the body siurface. 
Note that the body-surface values of c()q lie on the line 0 q = 17 in accordance with the 
zero-order boundary condition (eqs. (17)). However, the density function \j/Q and the 
longitudinal -velocity function Vq exhibit different types of singular behavior at the body 
surface for m 5^ 1; xJ/q goes to zero, and Vq goes toward minus infinity. This singu- 
larity at the body surface is an entropy -layer effect caused by the blunt nose of the body 
for m < 1. Van Dyke (ref. 16, page 186) notes the same singular behavior at the surface 
in a small -disturbance solution for hypersonic flow over a blunted wedge and observes 
that it occurs because the zero -order solution "is not a valid first approximation in the 
entropy layer." For the power-law body the effect is confined to a narrow region since 
the very high curvature in the nose area (infinite at x = 0) reduces the body slope 
rapidly. For example, when m = 0.80 and 6 = 0.5, the slope decreases from infinite 
to less than 1.0 before x/l = 0.006. 

Van Dyke (ref. 16) uses the method of matched asymptotic expansions to obtain a 
uniformly valid analytic solution for the blimted wedge. Adaptation of that method to the 
similarity -solution problem for power-law bodies probably would extend the solutions to 
the body surface; but, the application is complicated by the similar solutions being in 
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numerical form. Mirels (ref. 4) recognized the singularity at the body surface even 
though he did not calculate the longitudinal velocity function, which exhibits the singular 
behavior. He avoided the singularity by developing an approximate analytic solution and 
using it to obtain surface values of the zero-order functions (eqs. (20)). 

Shock-Displacement Constant 

O 

The variation of the calcvilated order -6 shock-wave displacement constant a 2 
with the power-law exponent m is shown in figure 5 for both two-dimensional and axi- 
symmetric bodies. The calculations made by using the iterative method for obtainii^ a 2 
gave essentially the same values as were obtained by the decomposition method (eq. (23)). 
However, the results from both methods are characterized by a singular discontinuity 
in a 2 which has no physical counterpart in the actual flow about power-law bodies. 

This singularity, which occurs near m = 0.817 for a = 0 and near m = 0.653 for 
a = 1 (two-dimensional and axially symmetric flow, respectively), represents a sudden 
decrease in the distance from the shock to the body as the power-law decreases, followed 
by a jump to a large distance across the discontinuity. Since this behavior is physically 
unrealistic, it must be produced by the mathematical processes (analytical and numerical) 
used to obtain the solutions. 

The upper part of figure 5 shows that the singularity is associated with a zero in 
the denominator of equation (23). Since there is a nearby zero in the numerator of this 
equation, it appears possible that the singularity in the quotient a 2 occurs because the 
zeros in the numerator and denominator are displaced relative to one another by accu- 
mulated errors in the ntunerical solution. (Coincidence of the zeros would make a 2 
mathematically indeterminate at that point but could allow a continuous variation of a 2 
with m, from which the value at the indeterminate point could be inferred.) To test this 
possibility, additional calculations were made by reducing the step size, reducing the 
range of the extrapolation, and, finally, extrapolating the whole of equation (23) rather 
than just the separate parts. The results of these calculations indicated that the position 
of the singularity and the values of a 2 for m greater than that at the singularity were 
virtually unaffected by the numerical procedure. The singularity was even unaffected by 
a complete reformulation of the problem in terms of momentum components instead of 
velocity components and by a change in extrapolation mode (ref. 14), although these 
changes did produce large changes in the computed value of a 2 for m less than that at 
the singularity (dashed line in fig. 5). Thus, removing the singularity would require a 
radical change in the mathematical process rather than simply changing parameters in 
the numerical integration and extrapolation schemes. For example, logarithmic terms 
might be required in the shock-wave and flow -variable expansions (eqs. (6) and (7)). 
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Shock displacement constant. 






The results of the present calculations should be reasonably good for body power- 

law exponents outside the range of significant influence from the singularity. Based on 

figure 5, the rai^e for good results is about 0.85 g m g 1.0 for the two-dimensional 

flow (a = 0) and 0.75 s rn ^ 1.0 for the axisymmetric flow (a = 1). Note that this is a 

more restricted range than that for which similarity solutions exist ( — - — < m ^ 1 

\3 + a 

(ref. 2)^) or for which entropy-layer effects are less than first order / — — < m ^ 1 
A \(3 + 0’)y + 2 

(ref. 17)) . 


Order-6 Functions 


9 

The order-6 similarity functions F 2 , 4^2’ ^2’ ^2 shown in figure 6 

for two-dimensional flow (a = 0) and in figure 7 for axisymmetric flow (a = 1). The 
results shown were obtained by using the decomposition technique for calculating the 
shock-displacement constant a 2 (eq. (23)), but essentially identical results were also 
obtained by using iteration. Numerical values of the functions are given in tables I and n. 
These order-6^ functions are seen to have some differences in behavior from the zero- 
order fimctions described previously. One obvious difference is that the curves describing 
these fimctions do not all emanate from a single point at r] = 1. The variation at r; = 1 
is due to the variation of the order-6^ boundary conditions at the shock with the body 
power law m and to the transfer of the boundary conditions from t] = ri„ to 77 = 1. 

o 

0 

All of the order-6'^ functions except the pressure function F« show singular 

^ 9 

behavior at the body surface for m < 1. These singularities in the order-6^ functions 
are a result of the entropy layer, just as in the case of the zero-order function In 

effect, the order-6 functions try to compensate for the large deviation of the zero-order 
longitudinal-velocity function at the surface, and so they become singular there also. The 
effect of the singularity is to make the order-6^ solution inapplicable at the body surface 
(except for the particular case of m = 1, for which the body does not have a blunt nose). 
However, the singularity should not affect the solution away from the surface, where the 
similarity functions are of order imity, so long as the constant a 2 is correctly deter- 
mined. Furthermore, the behavior of the order-6^ similarity function for the pres- 
sure Fo is quite regular all the way to the surface. Thus, the body surface pressure 

" O 

can be calculated to order-6 by using this function; however, the results must be sus- 
pect until checked against experiment or more exact results. 


It is because of the singular behavior of the similarity functions that the numerical 
integifttion cannot proceed all the way to the body surface. Neither do the extrapolations 
follow the singular functions in giving values at the surface; so, the calculated surface 
values of these fimctions do not represent the actual values of the singular functions. 
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which go to either plus or minus infinity at the surface. To the extent that they are useful 
for determining the value of a 2 , the extrapolations can be considered as providing a 
limiting process for this constant. 

Region of Validity of the Solutions 

Three basic assumptions were required in order to obtain the hypersonic similarity 
solutions for power-law bodies: (1) The body is slender enough that terms of order-6^ 
are negligible compared to unity; (2) the shock wave about the body is strong enough that 
terms of order e = are negligible compared to unity; and (3) the Mach number 

is large compared to tmity. (The second of these can be relaxed to « 1 if the first- 
order solution in e is applied.) However, even when these three assumptions are met 
overall for a particular power-law body, they generally are not all met in particular 
local regions. 

The first assumption, 6^ « 1, is obviously violated in the nose region of all blunt 
bodies, such as the power-law bodies for m < 1. Thus, the similarity solutions cannot 
be expected to apply at the nose of these bodies. But, the order-6^ solutions shoxold be 
particularly useful in providing an improved solution a moderate distance behind the nose. 
As mentioned in the previous section, the violation of the slender -body assumption by the 
blunt nose is also the cause of the singularities in the order-6^ functions at the body sur- 
face. Therefore, because the first assumption is violated at the nose, the order-6^ solu- 
tion does not apply at the body surface. 

The strong shock assumption, e « 1, is violated wherever the shock -wave angle 
approaches the Mach angle, sin“^(l/Moo). Unless the Mach number is extremely large, 
the shock wave will become weak far downstream from the nose df the body, and the simi- 
larity solutions will not apply in that region. 

The similarity solutions apply, then, in an intermediate region from behind the nose 
to somewhere in the vicinity of the base of the body, and, in the case of the order -6^ 
solution, only outside of the singularity at the body surface. The boundaries of this region 
depend on the Mach number and on the power-law exponent and slenderness of the body. 
Decreasing the Mach number or increasing the slenderness of the body tends to weaken 
the shock; decreasing the power-law exponent increases the nose bluntness but weakens 
the shock at the rear of the body. In any case, the boundaries of the region in which the 
solutions apply are not sharply defined but depend on the accuracy required in the results. 

Comparison With Other Solutions 

The only exact solutions available for comparison with the similarity solutions are 
those for flow over cones and wedges, corresponding to a power-law exponent of m = 1. 
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(A finite-difference calculation for m = 0.75 will be discussed in the next section.) 

Since the zero-order and order-6^ similarity solutions do not contain any Mach number 
dependence, the most appropriate comparison is at the hypersonic limit, — •». 

Results from the similarity solutions, in terms of the physical flow variables, are com- 
pared in figure 8 with the exact solutions for flow over a range of wedge angles at infi- 
nite Mach number. This flow, of course, is uniform behind the straight oblique shock 
wave, as indicated by the solutions for m = 1 in figures 3 and 6. In figure 9, results 
from the similarity solutions are compared with the exact solutions for flow over circular 
cones. Figures 8(a) and 9(a) show the variation of shock-wave angle with the body surface 
angle and include exact results from several Mach numbers in addition to M^ = «>. The 
other parts of figures 8 and 9 show the variations with the body surface angle of the pres- 
sure coefficient p/q^, the velocity V/u^ = and the velocity components 

u/u and v/u^ in the uniform flow behind the shock wave (fig. 8) and at the body sur- 
face (fig. 9). The similarity results are found from equations (6) and (7) with m = 1 as 
described in the appendix. The exact results are found from the oblique -shock relations 
in the two-dimensional case and from the charts of reference 15 and tables of reference 18 
in the axisymmetric case. 


Figures 8 and 9 show that for m = 1 the zero-order similarity solution agrees 
well with the exact solution for body surface angles up to about 0^ = 12°, whereas the 
order-6^ solution agrees well up to body aisles of about 0^ = 20°. As can be seen at the 
bottom of the figures, these cone or wedge angles correspond to slenderness-parameter 
values of about 6 =0.2 and 6 =0.4. The shock-wave angles calculated by similarity 
theory agree especially well with the exact solutions for infinite Mach number. In fact, 
over a substantial range of body angles the agreement is very good with the exact results 
for Mach numbers down to twelve. The similarity results for the magnitude of the veloc- 
ity also show good agreement for relatively large body angles (figs. 8(c) and 9(c)). Since 
the error is larger for the velocity components (figs. 8(d) and 9(d)), it must come mainly 
from error in predicting the direction of the velocity vector. This error in direction is 
shown in the upper part of figures 8(c) and 9(c), where it is compared to the curves 
6^0^ and -6^0]3- These curves represent the order of error expected from neglecting 
terms of order-6^ and of order-6'^, respectively. The error actually occurring is seen 
to be very close to that which was expected. It should be noted that this error in the 
direction of the velocity vector corresponds to an error in satisfying the boundary condi- 
tion that there is no flow through the body surface. That is, the velocity component nor- 


,-irb 


mal t£-the surface V sin tan - 0^1 should be zero. Figures 8(c) and 9(c) show 

^ o 9 

that this boundary condition is satisfied to order-6‘^ by the order-6 similarity solution 


for m = 1. 
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(b) Pressure p/q^o- 
Figure 8.- Continued 
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Figure 9.- Comparison of similarity solution for m = 1 with exact solution for 

flow over a circular cone at M = <». 
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(b) Surface pressure 

Figure 9.- Continued. 
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(d) Surface velocity components 
Figure 9.- Concluded. 



In figures 8(b) and 9(b), the pressure coefficients are compared also with the 
Newtonian prediction; p/q^o = 2 sin^e^. The Newtonian prediction is much more accu- 
rate for the conical flows (fig. 9(b)) than for the wedge flows (fig. 8(b)); however, even in 
the conical case the order-6^ similarity solution is closer to the exact solution for values 
of the similarity parameter 6 ^ 0.5. 

In figure 10, the variation of the flow variables from the shock to the body is shown 
for three values of the power-law exponent, m = 1.0, 0.85, and 0.75. The pressure, 
density, and velocity components, calculated from the axisymmetric zero-order and 
order-6^ similarity solutions, are shown at x/Z = 0.5 for a similarity parameter value 
of 6 = 0.4. Note that the zero-order results begin at = 1, which is the zero-order 
shock -wave location. The order-6^ results have a kink at r] = 1 which is caused by the 
switch from the Taylor series expansion to the integrated solution. 

The exact solution for a cone at = 20 with 0^ = 20° (ref. 18) is also shown 
for comparison in figure 10(a) with the case m = 1.0. Although the conical bodies are 
not exactly the same ^6^ = 20° corresponds to 6 = 0.3978), the nondimensionalized 
order-6^ similarity results agree well with the = 20 solution shown for all the 
variables except the density. As can be seen by the symbols representing the exact 
solution for at the shock wave, the density is the only one of the flow variables 

that is affected much by the differences between = 20 and — «>. The similarity 
solutions for the density agree exactly with infinite Mach number solution at the shock. 

The order -6^ solution for the other variables differs from the exact, infinite Mach number 
solution at the shock by amounts which are of order-6^, as expected from the approxima- 
tion to the oblique-shock relations used (eqs. (4)). On the other hand, the zero-order 

similarity solution is not accurate for a cone of this thickness; it is off by an amount of 
o 

order-6^, which is 16 percent for 6 = 0.4. 

For body power-law exponents other than m = 1 there is no exact solution avail- 
able for comparison with the order-6^ solution. However, there are some simple empir- 
ical methods for estimating the pressure on general bodies. Two of these methods will 
be used for comparisons. One is the Newtonian law Cp = 2 sin^e^^. As discussed by 
Hayes and Probstein (ref. 5), the Newtonian law corresponds to the limits y — 1.0 and 
M^ — °o; but, it is widely used for more general hypersonic flows in this or modified 
form. The other empirical prediction is the tangent-cone method, which takes as the 
pressure at any point on a body the pressure on the cone having the same surface angle 
as the body point. This method also is most accurate for M^ — since then the shock 
layer is very thin with little pressure change across it. Hayes and Probstein (ref. 5) 
give a thorough discussion of these two methods and their limitations. Only one limita- 
tion will be mentioned here; these methods give only the body surface pressure and are 
not complete flow-field solutions, as are the similarity solutions. 
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Exoct , Mco= 20. , 6[,= 20“ 


Zero-order similarity solutions 



(a) Power-law exponent, m = 1.0 (conical body). 

Figure 10.- Variation of flow variables from body to shock at x/f = 0.5 
for axisymmetric flow with 6 = 0.4. 
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(b) Power -law exponent, m = 0.85. 
Figure 10.- Continued. 
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(c) Power -law exponent, m = 0.75. 
Figure 10.- Concluded. 
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The calculated flow fields for bodies having m = 0.85 and m = 0.75 are shown 
in figures 10(b) and (c), respectively. In these cases the order-6^ solution is again a 
major correction to the zero-order solution. However, the singularities in some of the 
order-6^ similarity fimctions show up here at the body surface. Because of the singular- 
ities, the order-6^ values of the density and the velocity components are probably unreal- 
istic close to the surface. Fortunately, the pressure is well behaved all the way to the 
body surface, so that surface-pressure coefficients can be calculated. 

The surface pressures calculated from the similar solutions are seen in fig- 
ures 10(a), (b), and (c) to agree fairly well with the empirical predictions of the Newtonian 
and tangent-cone methods. However, as m goes from 1.0 to 0.85 to 0.75, the order-6^ 
similarity solution pressure drops faster than the tangent-cone and Newtonian pressures, 
so that the agreement becomes progressively worse. Whether the similarity solution or 
the empirical methods give a better representation of the actual pressure changes with 
body power law must be determined by comparison with experiment, as in the following 
section. 


Comparison With Experimental Results 

Only a limited amount of useful experimental data exist on the hypersonic flow fields 
about power-law bodies. These data consist mainly of measured shock-wave shapes and 
surface-pressure distributions for 3/4- and 2/3-power bodies of a few different fineness 
ratios. In this section the similarity -solution predictions for shock-wave shape and sur- 
face pressure for axisymmetric bodies with power-law exponent of m = 0.75 will be 
compared with the experimental results of Kubota (ref. 3) and Peckham (ref. 7). No 
comparison is made with results for smaller power-law exponents since valid order-6^ 
similarity solutions were not obtained in those cases. Also, no comparisons are made 
with experimental total -drag measurements (e.g., those of ref. 10) because of the uncer- 
tainty in calculating the skin-friction contribution. 

Shock shape . - Because the shock displacement constant a 2 is so small, the 
order-6^ shock-wave shape is only slightly different from the zero-order shape, an 
example of which is shown in figure 1 (m = 0.85; a = 0; 6 = 0.4). The zero-order and 

order-6 shock -wave shape predictions are compared to the shock-shape data for 
m = 0.75 from references 3 and 7 in figure 11. These data are presented in the corre- 
lation form discussed in the appendix to reference 7. The additional shock -wave -shape 
data from hypersonic flows over power-law bodies (refs. 8 and 9) are not presented here 
because they fall in ranges of very large x/D values, for which the order-6^ term of 
the shock -wave -shape equation is negligible. Also shown in figure 11 is the shock-wave 
shape, calculated by an asymptotic time -dependent finite -difference method (ref. 19, 
page 305), for hn 0.75-power body having a paraboloidal nose at = °o. 
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Figure 11.- Shock-wave shape correlation for axisymmetric 3/4-power bodies. 


The order-6^ similarity solution agrees very well with the finite-difference calcu- 
lation except at small x/D values where the latter is affected by its paraboloidal nose 
shape. The experimental data for the shock shapes are seen to correlate well with one 
another but to fall somewhat above the similarity solution and finite -difference predic- 
tions. This disagreement between experiment and theory is due in part to the disagree- 
ment in Mach numbers (M^ = 6.85 and 7.7 for experiment, Moo = °° for theory). It is 
comparable to the shift in shock location with Mach number for cones (fig. 9(a)). An 
additional cause for the disagreement between the experiment and theory is the outward 
displacement of the flow by the growth of the viscous boundary layer on the experimental 
bodies. 

The effect of the order-6^ term in the similarity solution for the shock -wave shape 
is seen in figure 11 to increase as x/G decreases. This is expected since small values 
of x/D correspond either to small values of x/l or to large values of D/f, that is, 
either to points near the nose of the body (where the slope is large) or to bodies which 
are less slender (and thus have larger 6-values). Note that the agreement of experiment 
with theory is better in this region of smaller x/D values, as would be expected since 
the shock location is closer to that for infinite Mach number for larger body slopes. The 
slope of the order-6^ similarity solution agrees very well with the slope of Kubota's data 
for D/r = 0.0555 (see circles in fig. 11); this agreement in the slope on a log-log plot 
indicates good agreement with the power-law exponent of the physical shock-wave shape. 

Pressure distribution .- The experimental pressure distributions obtained by 
Kubota (ref. 3) and Peckham (ref. 7) are shown in figure 12 for the same power-law 
bodies (m = 0.75) as used for the shock-wave shapes. For comparison, the zero-order 
and order-6^ similarity -solution predictions are also shown, along with the two empirical 
pressure distributions. For the three bodies having fineness ratios i = lj2r^{X) of 
about 2, the similarity solutions, as well as the two empirical methods, give pressure 
distributions in good agreement with the experimental data. 

2 

For these cases the order-6 similarity solution is very nearly the same as the 
zero-order solution except at the front of the body. The Newtonian prediction falls 
slightly higher than the similarity solutions behind the nose region, but curves repre- 
senting these three methods are below the data points. Since at hypersonic speeds the 
viscous boundary layer tends to displace the flow outward, raising the pressure above 
that which would occur for inviscid flow, the theoretical inviscid pressure levels are 
ejq)ected to fall slightly below those actually measured. For example, by applying a 
boundary -layer displacement correction to the zero-order similarity solution for the 
pressure on his fineness -ratio -2. 13 body, Kubota (ref. 3) obtained excellent agreement 
with his experimental data. (See squares in fig. 12.) Since the tangent-cone pressure 
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distribution falls slightly above the experimental data for the three higher fineness-ratio 
bodies, the other methods are somewhat preferable. 

Only in the case of the body having a fineness ratio of nearly 1, however, does a 
real difference between the methods appear in figure 12. In particular, the difference 
between the zero-order and order-6^ pressure distributions becomes substantial. The 
zero-order result lies above the experimental data by an amount which increases rapidly 
toward the front of the body. On the other hand, the order-6^ result lies below the data 
by an amount which, for most of the body, is only moderately larger than that for the 
finer bodies. This amount is on the order of the expected boundary -layer displacement 
effect. It is only at the front of the body that the order-6^ result begins to diverge 
markedly from the experimental pressure distribution. Since the value of the slender- 
ness parameter is 5 = 0.555 for this case, it is not surprising that the order-6^ simi- 
larity solution should begin to fail as the body surface angle increases at the front of the 
body. (Kubota's value 6 = 0.485 shown with his data corresponds to 617^ = 0.485 as 
used herein.) This is about the same value of 6 as the limit for good results in the 
wedge and cone cases (figs. 8 and 9). The Newtonian method gives excellent agreement 
with the experimental data in this case, but this must be somewhat fortuitous in that no 
correction was made to account for the boundary -layer displacement or body -curvature 
effects. The tangent-cone method again lies somewhat above the experimental data. 

CONCLUSIONS 

By beginning with the equations for conservation of mass, conservation of momentum, 
and conservation of energy for the inviscid, two-dimensional or axisymmetric adiabatic 
flow of an ideal gas, similarity solutions have been found which give the flow field to 
order-6^ about power-law bodies in the hypersonic limit M^o — °°, where 6 is a body 
slenderness parameter and Moo is the free-stream Mach number. On the basis of this 
investigation the following conclusions can be made: 

1. Order-6^ solutions were obtained which are independent of the slenderness 
parameter 6. Thus, the functions expressing the solutions are universal in that they 
apply for all values of 6 for which 6'^ « 1. The relations between those similarity 
functions and the physical flow variables are relatively simple. 

2. In the present formulations the value of ag, the shock displacement constant in 
the order-6^ solution, goes through plus and minus infinity at about m = 0.817 in the 
two-dimensional case and about m = 0.653 in the axisymmetric case, where m is the 
body power-law exponent. Because the singularity does not correspond to actual flow 
conditions, it must arise through the mathematical development. Since the singularity 
was not removed by any of the variations in solution procedure tried, the present results 
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are limited to a range judged relatively free of effects from the singularity 

(0.85 g m ^ 1.0 for two-dimensional flow; 0.75 g m g 1.0 for axisymmetric flow). 

3. In comparisons with the exact solutions for inviscid flow over wedges and circu- 
lar cones, the order-6^ similarity results give excellent agreement for 6 less than 
about 0.4, corresponding to we(%e or cone angles up to about 20°. Over an even larger 
range, the order-6^ surface -pres sure predictions were superior to the Newtonian pres- 
sure law. The order-6^ results were a significant improvement over the zero-order 
results for body angles greater than about 12°. 

4. In comparisons with experimental shock-wave shapes and surface -pressure dis- 
tributions for 3/4 -power axisymmetric bodies, the order -6^ similarity solutions gave 
good results, considering that Mach number and boundary -layer displacement effects are 
not included in the theory. For body fineness ratios near 2, the effects of the order-6^ 
terms are significant only very near the body nose; whereas for a fineness ratio near 1, 
the order-6^ terms had a large effect over almost the entire body. These good results 
for the surface pressure were obtained despite the singular behavior of some other vari- 
ables at the surface. 

5. Although the order-6^ similarity solutions were developed for the hypersonic 
limit — °o, the derivation shows that they are compatible with order-e solutions, 
where e = 1^(M^6)^, which introduce Mach nvunber effects. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
July 1, 1975 
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APPENDIX 


APPUCATION OF THE ORDER- 6^ SOLUTIONS 

In order to facilitate application of the order-6^ similarity solutions, the equations 
relating the similarity functions to the physical -flow variables are presented here. It is 
assumed that the body fineness ratio f and power-law exponent m are known, so that 
the body shape is 

% ^ 

I 2f\^f ) 

The value of cr is 0 for a two-dimensional body and is 1 for an axisymmetric body. 
Then, the values of 77^, the body-surface similarity coordinate, and of a^, the shock 
displacement constant, are found in table III for the given m. The value of 6, the 
slenderness parameter, is calculated from 

6 = 

2f77b 

The order-6^ shock-wave shape is 


R 

r 




-2(l-m) 


The similarity coordinates of any point (x,f) in the flow field between the body and the 
shock wave are 


By taking the values of the similarity functions from table I (for two-dimensional flow) 
or table II (for axisymmetric flow), the pressure, density, and velocity components at 

o 

any point are given to order -6 by 

J-= 26^|FQ(77)m^|"^^^““) + 6^F2(T7)m'^^"‘^^^"'"^ 
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P OO 

i = 1 + 6^„(n)m2{-2(l-“) 
u ^ 

OO 


Because the entropy layer makes the similarity solutions for the velocity and density 
inapplicable at the body surface for m < 1, these flow variables must be found by other 
means. The similarity solution for the pressure, on the other hand, does give a reason- 
able result at the surface, which can be used to help in calculating the density and 
velocity. The basis of the method for calculating the density and velocity at the surface 
is the fact that the streamline along the surface must have passed through the normal 
shock wave at the blunt nose of the body. After passing through the shock, the flow is 
isentropic along streamlines. Thus, by letting the subscript n denote conditions just 
behind the normal shock, and by using isentropic and normal shock relations (e.g., ref. 14), 
the density at the surface is 


Poo PcoPn Poo\Pn/ 



1/y 


' (r + 1)mJ " 

(y - 1)M^2 + 2 


y(y + 1)mJ^ 
2yUj - (y - 1) 


-ii/y 




Similarly, the magnitude of the velocity at the surface is given by 
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or 



2y 


y - 1 




2y P°o 

Tt 

y - 1 p 

^ OO 

Too 







In the limit as the density and velocity become 


and 





1/y 



Here 


Pb 


_ P(^>^b) 

6^p u 2 

^ oo oo 


Since the flow at the 


is the normalized surface pressure given by the similarity solutions, 
body is parallel to the surface, the velocity components are 
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TABLE I.- SIMILARITY FUNCTIONS FOR TWO-DIMENSIONAL FLOW (cr = 0) 



— 



- • ■ 

— 


1 


^ 1 

^0 

^0 

''o 

^0 

^2 

^2 

'’2 j 

♦2 





m = 1.00 





1 *00000 

•83333 

6*0000 

.83333 

•83333 

-.5556 

0.0000 

-.5556 

-.6944 

•99219 

•83333 

6.0000 

•83333 

.83333 

-•5556 

0.0000 

-.5556 

-.6944 

•98438 

•83333 

6*0000 

.83333 

•83333 

-.5556 

0.0000 

-.5556 

-.6944 

.97656 

•83333 

6*0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

.96875 

.83333 

6.0000 

•83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•96094 

•83333 

6.0000 

•83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•95313 

.83333 

6*0000 

.83333 

.83333 

-.5556 

0.0000 

-•5556 

-.6944 

•94531 

.83333 

6*0000 

• 83333 

•83333 

-.5556 

0.0000 

-.5556 

- .6944 

•93750 

.83333 

6*0000 

.83333 

.63333 

-.5556 

0.0000 

-.5556 

-.6944 

•92969 

• 83333 

6*0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•92188 

• 83333 

6*0000 

•83333 

•83333 

-.5556 

0.0000 

-.5556 

-.6944 

•91406 

.83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•90625 

.83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

- .6944 

•89844 

.83333 

6.0000 

•83333 

•83333 

-.5556 

0.0000 

-.5556 

-.6944 

•89063 

•83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•88281 

•83333 

6.0000 

.83333 

.83333 

-.5556 

0*0000 

-.5556 

-.6944 

•87500 

•83333 

6 . 0 O 00 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

• 86719 

•83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

•85938 

.83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

.85156 

.83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

.84375 

.83333 

6.0000 

.83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 

.63594 

.83333 

6.0000 

•63333 

•83333 

-.5556 

0.0000 

-.5556 

-•6944 

.83333 

•63333 

6.0000 

•83333 

.83333 

-.5556 

0.0000 

-.5556 

-.6944 





m = 0.95 





1 .00000 

•83333 

6.0000 

.83333 

•83333 

-.6414 

-.8290 

-.6107 

-.7259 

•99219 

•83030 

5.9629 

.83373 

•83232 

-.6307 

-.7645 

-.6086 

-.7232 

•98438 

.82735 

5.9254 

.83425 

.83133 

-.6205 

-.7018 

-.6060 

-.7207 

•97656 

•82448 

5.8876 

•83489 

•83038 

-.6107 

-.6409 

-.6054 

-.7185 

•96875 

•82170 

5.8493 

.83566 

•82944 

-.6014 

-.5813 

-.6042 

-.7165 

•96094 

•81899 

5.8103 

.83658 

• 82854 ^ 

-.5924 

-.5231 

-.6034 

-.7140 

•95313 

•81636 

5.7706 

.83764 

•82765 

-.5839 

-.4658 

-.6029 

-.7132 

•94531 

•81380 

5.7300 

•83687 

.82679 

-.5757 

-.4094 

-.6028 

-.71 19 

.93750 

.81 131 

5.6864 

•84029 

.82595 

-.5678 

-.3536 

-.6031 

-.7109 

.92969 

•80888 

5.6454 

•84191 

.82513 

-.5602 

-.2982 

-.6038 

-.7101 

•92188 

•80652 

5.6010 

•84375 

.82434 

-.5529 

-.2430 

-.6049 

-.7095 

•91406 

•80422 

5.5548 

•84586 

•82356 

-.5460 

-.1875 

-.6065 

-•7093 

.90625 

.80199 

5.5064 

.84826 

•82280 

-.5392 

-.1316 

-.6087 

-.7094 

.89844 

•79981 

5.4554 

•85101 

•82205 

-.5328 

-.0749 

-.61 15 

-.7098 

.89063 

.79769 

5.4013 

•85418 

.82133 

-.5266 

-.0167 

-.6151 

-.71 06 

•86281 

•79562 

5.3433 

•85786 

•82062 

-.5206 

• 0434 

-.6196 

-.71 19 

•87500 

•79361 

5*2805 

•86216 

.81993 

-.5148 

• 1062 

-.6251 

-.7138 

.86719 

.79166 

5.21 1 3 

•86727 

.81925 

-•5093 

• 1730 

-•6322 

-.7164 

.85938 

•78977 

5 '. 1 338 

•87346 

•81858 

^.5039 

• 2453 

-.641 1 

-.7199 

.85156 

•78793 

5.0448 

•88116 

•81 793 

-.4988 

.3257 

-.6528 

-.7247 

.84375 

.78616 

4.9387 

•89113 

•81730 

-•4939 

■ 4185 

-•6686 

-.7314 

.83594 

.78445 

4*8049 

•90489 

•81668 

-.4892 

• 5324 

-.6916 

-.7412 

•82813 

• 78282 

4*6180 

.92629 

•81607 

-.4847 

• 6879 

-•7293 

-.7573 

•62031 

•781 30 

4.2817 

• 97 M 6 

•81547 

-.4805 

.9657 

-.8157 

-•7925 

•81507 

.78040 



•81507 

-.4780 
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TABLE I.- Concluded 


n 

L_jo 

T ' , 

1 ^0 

''0 


1 .. 

^2 

1 

♦2 

1 .00000 

.83333 

6.0000 

] 

•83333 

m = 0.90 

•83333 

-.7510 

-1 .7618 

-.6865 

— • 7681 

#992 1 9 

•82693 

5.9221 

.83417 

•831 19 

-.7281 

- 1 .6059 

-.6836 

-*7624 

•98436 

#82075 

5.8444 

•83523 

.82910 

7057 

- 1 .4558 

-.6814 

-.7572 

#97656 

•81478 

5.7668 

.83652 

.82707 

-.6844 

- 1.3112 

-.6799 

-*7525 

.96875 

.80900 

5.6893 

•83807 

•82508 

-.6642 

- 1 .1716 

-.6792 

-.7481 

•96094 

#80342 

5.61 1 5 

.83989 

.82315 

-.6449 

-1 .0365 

-.6793 

-.7442 

#95313 

#79802 

5.5335 

•84200 

.82126 

-.6266 

-.9056 

- . 68 02 

-* 740 e 

.94531 

#79261 

5.4548 

.84441 

.81941 

-.6091 

-.7783 

-.6820 

-.7378 

•93750 

#78776 

5.3755 

.84716 

.81761 

-.5926 

-•6543 

-.6846 

-.7352 

#92969 

. 78288 

5.2952 

.85028 

.81585 

-.5767 

-.5331 

-.6883 

-.7331 

#92188 

#77816 

5.2137 

•85380 

•81414 

-.5617 

-.4143 

-.6931 

-.7315 

.91406 

#77360 

5. 1 307 

.85777 

•61246 

-.5473 

-.2974 

-.6991 

-*7304 

.90625 

.76919 

5.0458 

.86225 

♦ 81 oai 

-.5337 

1820 

-.7064 

-•7299 

.89844 

#76493 

4.9503 

.86731 

.80921 

-.5206 

-.0675 

-.7153 

-.7300 

.89063 

#76082 

4.8691 

.87303 

.80764 

-.5082 

.0466 

-.7260 

-.7307 

.88281 

•75685 

4,7761 

.87953 

•80610 

-.4964 

.1612 

-.7388 

-*7322 

.87500 

• 75302 

4.6792 

.88693 

.80459 

-.4851 

.2768 

-.7541 

-.7345 

.86719 

.74933 

4 . 5776 

.89543 

*8031 1 

-.4744 

• 3946 

-.7725 

-.7377 

.85938 

. 74578 

4 .4700 

.90527 

.00167 

-.4641 

.5157 

-.7948 

-*7422 

.85156 

• 74237 

4.3551 

.91681 

-80025 

-.4544 

.6418 

-.8220 

-.7481 

.84375 

.73910 

4.2308 

.93053 

.79886 

-. 4452 

.7750 

-.8559 

-.7558 

. 8359<4 

• 73598 

4.0941 

.94718 

.79749 

-.4365 

.9104 

-.8991 

-.7656 

.82613 

• 73301 

3.9406 

.96798 

.79615 

-.4282 

1 .0769 

-.9560 

-.7793 

.82031 

• 73020 

3.7631 

.99502 

.79483 

-.4205 

1 .2587 

-1 .0347 

-.7977 

•81250 

• 72757 

3.5483 

1 .03251 

.79353 

-.41 33 

1 .4788 

- 1 .1526 

-.8244 

.80469 

.72515 

3.2672 

I .09065 

.79226 

-.4066 

1 .7720 

-1 .3558 

-.8674 

#79688 

• 72301 

2.8296 

1.20714 

.79100 

-.4007 

2.2552 

- 1 .8406 

-.9565 

.78988 

.72154 



• 78988 

-.3966 





m = 0.85 


1 .00000 

.83333 

6.0000 

.83333 

.83333 

-.8590 

- 1 .6653 

-.7974 

-.8331 

•99219 

•82320 

5.6770 

.83464 

•82993 

-.8221 

- 1.3962 

-.7965 

-.8145 

•96436 

.81345 

5 . 7 S 59 

.83627 

•82660 

-.7872 

- i . 1 400 

-#7968 

-.8066 

.97656 

.80409 

5.6367 

.83823 

.82334 

-.7544 

- .8955 

-.7984 

-.7992 

.96875 

• 79509 

5.5190 

.84055 

.82016 

-.7234 

-.6620 

-.8014 

-.7925 

•96094 

• 78644 

5.4028 

.84325 

.81 704 

-.6942 

-.4386 

-.8058 

-•7863 

.95313 

.77813 

5.2879 

.84634 

.81400 

-.6667 

-•2244 

-. 81 1 7 

-#7808 

.94531 

• 770 1 4 

5 . 1 739 

.84986 

.81 1 02 

-.6407 

-.0107 

-. 81 92 

-.7758 

.93750 

* 76246 

5.0609 

.85383 

.8081 1 

-.6162 

. 1 793 

-.0285 

-.7715 

.92969 

• 75509 

4 .9484 

•85829 

.80526 

-.5932 

.3703 

-.8397 

-#7678 

*92188 

* 74800 

4.8364 

.86328 

.80247 

-.5714 

• 5552 

- .8530 

-#7647 

.91406 

.741 19 

4.7247 

•86886 

.79974 

-.5509 

.7345 

— # 8666 

-.7623 

.90625 

* 73466 

4.6128 

.87507 

•79706 

-.5316 

.9092 

-.8867 

-.7606 

.89844 

.72838 

4.5007 

•88199 

♦79444 

-.5134 

1 .0799 

-.9078 

-.7597 

.89063 

.72237 

4.3880 

•88970 

.79187 

-.4963 

1 .2473 

-.9322 

-.7595 

.88261 ; 

.71660 

1 4.2745 

1 •89629 

.78935 

480 1 

1.4124 

-.9605 

-.7603 

.87500 

.71107 

4. 1596 

i .90789 

•70687 

-.4650 

! 1 .5759 i 

-#9932 

-.7619 

.06719 

1 *70577 

1 4.0432 

•91864 

.78445 

-.4507 

1 .7389 

- 1 .031 1 

-.7645 

.85938 

.70071 

3.9246 

.93071 

.78207 

-.4373 

1 .9023 

- 1 .0754 

-#7684 

*85156 

*69588 

3.8034 

.94435 

•77973 

1 -.4248 

2.0674 

- 1# 1274 

-.7735 

• 8^(375 

.69127 

3%6759 

.95984 

.77743 

-.4131 

2.2355 

- 1 • 1889 

-.7802 

.83594 

• 68688 

3.5504 

.97757 

.77517 

-.4021 

2.4085 

-1 .2624 

-.7886 

.82813 

.68271 

3.4160 

•99805 

.77295 

-.3919 

2.5885 

- 1 .3514 

-.7993 

.82031 

* 67877 

3.2770 

1 .02200 

.77076 

-.3825 

2.7781 

- 1.4610 

-.8126 

.81 250 

* 67504 

3. 1292 

1 .05040 

.76861 

-.3737 

2.9612 

- 1 .5988 

-#8294 

.80469 

*67155 

2.9714 

1 .08478 

•76649 

-.3658 

3.2029 

- 1 .7770 

-.8507 

*79688 

. 66829 

2.8002 

1 % 12745 

. 76440 

-.3585 

3.4510 

- 2.0159 

-#8763 

.78906 

.66527 

2.61 1 0 

!• 18236 

•76233 

-.3520 

3.7377 

- 2# 3529 

-.9151 

.78125 

* 66252 

2 .3959 

1 .25684 

.76029 

-.3463 

4.0845 

- 2.8652 

— # 9666 

• 77344 

• 66006 

2 . 1 407 

!• 36668 

.75827 

-.3416 

4.5347 

- 3*7427 

- 1 .0446 

.76563 

* 65794 

1.8131 

1 .55607 

•75628 

-.3379 

5. 1985 

- 5.6225 

- 1 . 1824 

.75781 

%65627 

1 .3010 

2.04995 

•75429 

-.3358 

6*5593 

- 13.0166 

- 1 .5489 

.7531 1 

.65566 



.7531 1 

-.3364 


1 



TABLE n.- SIMILARITY FUNCTIONS FOR AXISYMMETRIC FLOW (ff = 1) 


n 

1 Fq 

'*'0 


*0 

^2 

^2 

^2 _ 

*2 




m = 1.00 


1 .00000 

.83333 

6.0000 

.83333 

.83333 

-.6091 

.4248 

_ 

.6091 

-.6799 

.99219 

.84055 

6.0371 

.84085 

.84088 

-.6306 

.3445 

- 

.6250 

-.6959 

.98438 

.84705 

6.0704 

.84823 

.84835 

-.6495 

.2746 

- 

.6401 

-.71 12 

.97656 

.85284 

6.1000 

.85550 

.85576 

-.6660 

.2142 

- 

.6546 

-.7259 

.96875 

.85796 

6.1261 

.86269 

.86315 

-.6804 

.1622 

- 

.6685 

-.7402 

.96094 

.86240 

6.1487 

.86981 

.87053 

-.6927 

.1180 

- 

.682 0 

-.7542 

.95313 

. 866 1 7 

6.1679 

.87689 

.87793 

-.7030 

.0812 

- 

.6951 

-.7679 

.94531 

.86927 

6.1837 

.88394 

.88536 

-.7115 

.0513 

- 

.7080 

-.7814 

.93750 

.87169 

6 . 1 960 

.89098 

.89284 

-.7180 

.0283 

- 

.7206 

-.7949 

.92969 

.87342 

6.2048 

.89804 

•90040 

-.7226 

.0120 

- 

.7332 

-.8084 

.92188 

.87445 

6.2100 

.90512 

.90804 

-.7254 

.0025 

- 

.7457 

-.8218 

.91494 

.87487 

6.2131 

.91142 

.91488 

-.7270 

-.0032 

- 

.7567 

-.8338 




m = 0.95 




1 .00000 

.83333 

6.0000 

.83333 

.83333 

-.6542 

.0228 


.6390 

- . 6959 

. 992 1 9 

.83754 

5.9994 

.84131 

.83988 

-.6651 

.0013 

- 

.6530 

-.7093 

.98438 

.84113 

5.9935 

.84939 

.84640 

-.6742 

-.0103 

- 

.6669 

-.7224 

.97656 

.84415 

5.9819 

.85765 

.85292 

-.6815 

-.0125 

- 

.6808 

-.7354 

.96875 

.84659 

5.9640 

. 866 1 5 

.85945 

-.6870 

-.0054 

- 

.6949 

-.7485 

.96094 

.84846 

5.9392 

.87500 

.86601 

-.6911 

.0115 

- 

.7095 

-.7618 

.95313 

.84978 

5.9060 

•88433 

.87262 

-.6935 

.0388 

- 

.7248 

-.7755 

.94531 

.85053 

5.8624 

.89437 

.87929 

-.6945 

.0781 

- 

.7413 

-.7898 

.93750 

.85073 

5.8047 

.90549 

•88604 

-.694 1 

.1324 

- 

.7597 

-.8051 

.92969 

.85036 

5.7255 

.91844 

.89288 

-.6922 

.2084 

- 

.7815 

-.8222 

.92188 

.84943 

5.6056 

.93522 

.89983 

-.6889 

.3232 

- 

.81 05 

-.8429 

.91406 

•84796 

5*3556 

.96541 

.90692 

-.6843 

.5581 

- 

.8658 

-.8760 

.91 034 

.8471 1 



.91034 

-.6817 








m = 0.90 






1 .00000 

.83333 

6.0000 

.83333 

.83333 

-.7075 

-.4384 

_ 

.6750 

-.7156 

.99219 

.83420 

5.9581 

.84181 

.83877 

-. 7063 

-.3849 

- 

.6876 

-.7260 

.98438 

.83462 

5.9101 

.85064 

.84424 

-.7040 

-.3240 

- 

.7009 

-.7366 

.97656 

.83462 

5.8554 

.85992 

.84974 

-.7008 

-.2553 

- 

.7151 

-.74 77 

.96875 

.83421 

5.7931 

.86978 

.85531 

-.6966 

-.1783 

- 

.7303 

-.7593 

.96094 

.83338 

5.7220 

.88036 

.86094 

-.6914 

-.0921 

- 

.7471 

-.7716 

.95313 

.83215 

5.6403 

.89192 

• S 6666 

-.6854 

.0050 

- 

.7660 

-.7849 

.94531 

.83051 

5.5452 

.90479 

.87246 

-.6786 

• 1 154 

- 

.7877 

- . 7995 

.93750 

.82848 

5.4321 

.91957 

.87837 

-.6709 

• 2434 

- 

.81 37 

-.8159 

.92969 

.82606 

5.2928 

.93733 

.88440 

-.6624 

.3970 

- 

.8464 

-.8353 

.92188 

.82328 

5.1101 

.96043 

.89056 

-.6532 

.5932 

- 

.8918 

-.8597 

.91406 

.82018 

4.8352 

.99587 

• 896 B 6 

-.6434 

.8813 

- 

.9676 

-.8956 

.90625 

.81691 

4.1 357 

1 .09932 

.90331 

-.6334 

1*61 

-1 

.2321 

-.9935 

.90465 

•81630 



•90465 

-.6315 
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TABLE II.- Concluded 



( ^0 


I ^0 
m = 0.85 


.83333 

.83048 

.82743 

.82415 

.82067 

.81698 

.81309 

.80901 

.80474 

.80030 

.79572 

.79105 

.78640 

.78214 

.78174 


6.0000 

5.9124 

5.8191 

5.7195 

5.6125 

5.4968 

5.3705 

5.2312 

5.0746 

4.8944 

4.6790 

4.4035 

3.9979 

2.9087 


.83333 
.84235 
.85198 
.86233 
.87357 
.88590 
.89961 
.91515 
.93319 
.95489 
.98241 
1 .02070 
I .08536 
1.34315 


.83333 

.83752 

.84180 

.84617 

.85064 

.85521 

.85991 

.86472 

.86967 

.87476 

.87999 

.88539 

.89095 

.89668 

.89743 

m = 0.80 


-.7700 

-.7546 

-.7392 

-.7239 

-.7086 

-.6934 

-.6782 

-.6630 

-.6478 

-.6326 

-.6176 

-.6027 

-.5883 

-.5754 

-.5742 


-.9473 
-.7969 
-.6442 
-•4882 
-.3278 
-.1615 
.0128 
.1980 
. 3984 . 
.6214 
.8802 
1 .2035 
1 .6771 
3.0601 


-.7188 
-.7310 
-.7447 
-.7603 
-.7782 
-.7991 
-.8236 
-.8531 
-,8896 
-.9367 
- 1.0014 
- 1 ,1007 
-1 .2928 
- 2,3746 


.83333 

6.0000 

.83333 

.83333 

-.8390 

- 1 .4237 

-.7712 

.82633 

5.8617 

.84295 

.83612 

-.8067 

- 1.1475 

-,7842 

.81941 

5.7196 

.85341 

.83904 

-.7762 

-.8785 

-.8000 

.81256 

5.5731 

.86487 

.8421 1 

-. 7472 

-.6149 

-.8190 

.80578 

5.4212 

.87751 

.84533 

-.7197 

-.3547 

-.8419 

.79906 

5.2627 

.89156 

.84869 

-.6934 

-.0957 

-.8695 

.79239 

5.0960 

.90735 

.85220 

-.6683 

.1647 

-.9029 

.78579 

4.9192 

.92535 

.85587 

-.6443 

.4295 

-.9441 

.77925 

4.7292 

.94626 

.85970 

-.6214 

.7030 

-.9957 

.77279 

4.5216 

.97114 

.86370 

-.5996 

.9913 

- 1 .0623 

.76644 

4.2904 

1.00183 

.86786 

-.5788 

1 .3041 

- 1 . 1520 

.76023 

4.0232 

1.04173 

.87219 

-.5591 

1 .6582 

- 1.2812 

.75426 

3.6971 

1 .09836 

.87671 

-.5407 

2.0890 

- 1 .4890 

.74865 

3.2526 

1. 19402 

.88140 

-.5240 

2.6922 

- 1 .9067 

.74380 

2.3985 

1 .47814 

.88628 

-.5102 

4.0109 

- 3.6855 

.74265 



.88798 

-.5073 




m = 0.75 


.83333 
•82165 
.81044 
•79967 
•78932 
.77937 
.76979 
.76058 
•75172 
• 74322 
. 73509 
.72735 
.72003 
.71322 
.70702 
.70168 
.69806 


6.0000 

5.8050 

5.6103 

5.4150 

5.2183 

5.0190 

4.8161 

4.6082 

4.3932 

4.1688 

3.9316 

3.6765 

3.3956 

3.0747 

2.6846 

2.1423 


.83333 

.83333 

-.8898 

- 1.4529 

•84360 

.83452 

-.8391 

- 1.0224 

.85495 

.83590 

-.7922 

-.6147 

.86755 

.83747 

-. 7487 

-.2265 

•881 57 

.83923 

-.7085 

• 1456 

•89729 

.84119 

-.6711 

.5048 

.91502 

.84333 

-.6365 

.8546 

.93524 

.84566 

-.6044 

1 . 1 989 

.95858 

.84817 

-.5746 

1 .5421 

.98598 

.85088 

-.5470 

1 .8896 

1 .01886 

-85377 

-.5217 

2.2491 

1 .05952 

•85686 

-.4984 

2.6315 

1 . 1 1208 

.86013 

-.4773 

3.0546 

1 . 18477 

.86360 

-.4586 

3.5515 

1 .29784 

.86726 

-.4426 

4.1965 

1 .52579 

.871 10 

-.4301 

5.2307 


.87507 

-.4253 







-.8253 
-.8409 
-.8609 
-.8857 
-.9162 
-.9536 
-.9997 
- 1 .0566 
- 1 .1280 
-1 .2193 
-1 .3391 
-1 .5030 
- 1#7407 
- 2.1186 
- 2.8255 
- 4.7293 


TABLE m.- SHOCK DISPLACEMENT CONSTANT AND 
PRESSURE FUNCTIONS AT THE BODY SURFACE 


m 

% 

a = 0 

F^Cnfa) 

^2 

1 .00 

.83333 

.83333 

-.55556 

• 16667 

.95 

.81507 

.78040 

-.47804 

• 17500 

.90 

.78988 

•72154 

-.39659 

.17618 

.85 

.7531 1 

.65566 

a = 1 

-. 3364 1 

. 10485 

1 .00 

.91494 

.87487 

-.72701 

.08496 

.95 

.91034 

.8471 1 

-.68169 

.0866 1 

.90 

.90465 

.81630 

-.63154 

,08768 

.85 

.89743 

.78174 

-.57424 

,06705 

.80 

.88798 

. 74265 

— • 50734 

.08135 

.75 

.87507 

.69806 

-.42532 

.05812 
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